Interplay of superradiance and disorder in the Anderson Model 
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Using a non-Hermitian Hamiltonian approach to open systems, we study the interplay of disorder 
and superradiance in a one-dimensional Anderson model. Analyzing the complex eigenvalues of the 
non-Hermitian Hamiltonian, a transition to a superradiant regime is shown to occur. As an effect of 
openness the structure of eigenstates undergoes a strong change in the superradiant regime: we show 
that the sensitivity to disorder of the superradiant and the subradiant subspaces is very different; 
superradiant states remain delocalized as disorder increases, while subradiant states are sensitive to 
the degree of disorder. 
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I. INTRODUCTION 

Nanoscopic systems in the quantum coherent regime 
are at the center of many research fields in physics, rang- 
ing from quantum computing and cold atoms to trans- 
port in nanoscale and mesoscopic systems. Transport in 
the quantum coherent regime can be considered one of 
the central subjects in modern solid state physics [TJ [5] 
and in cold atom physics [3]. Transport properties de- 
pend strongly on the degree of openness of the system. 

As a consequence of quantum coherence many inter- 
esting features arise. Here we focus on two important 
effects induced by quantum coherence: Anderson local- 
ization [4] and Dicke superradiance [5]. Anderson Lo- 
calization is driven by intrinsic disorder and consists in a 
suppression of diffusion due to an exponential localization 
of the eigenfunctions of the system. Dicke superradiance 
is driven by the fact that the system is open, namely 
coupled to an external environment characterized by a 
continuum of states. To explain the superradiance ef- 
fect, consider a discrete quantum system coupled to an 
environment having a continuum of states. The system- 
environment coupling alters the unperturbed energy lev- 
els: it causes an energy shift and the appearance of a res- 
onance width (inverse lifetime) for each level. For weak 
coupling strength, all resonance widths are roughly equal. 
However, once the coupling strength reaches a critical 
value, the widths start to overlap, and width segregation 
occurs. In this regime, almost the entire (summed up) 
decay width is allocated to just a few short-lived "super- 
radiant" states, while all other states are long-lived (and 
effectively decoupled from the environment). We call this 
segregation the "superradiance transition" (ST). 

In the superradiant regime, the effect of the opening 
is large, and cannot be treated perturbatively. Thus, a 
consistent way to take the effect of the opening into ac- 
count for arbitrary coupling strength between the system 
and the outside world is highly desirable. The effective 
non-Hermitian Hamiltonian approach to open quantum 



systems has been shown to be a very effective tool in ad- 
dressing this issue [5J . Non-Hermitian Hamiltonians have 
been already employed to study realistic open quantum 
systems, such as transport through quantum dots [7], 
superradiance in cold atoms [3j and nuclear physics [8]. 
The superradiance effect has also been studied using ran- 
dom matrix theory [pj [ID], in microwave billiards [TT] 
and in paradigmatic models of coherent quantum trans- 
port P21 [13]. As an example of the importance of the 
ST, maximum transmission in a realistic model for quan- 
tum transport was shown to be achieved exactly at the 

st [ig. 

In this paper we analyze a one-dimensional Anderson 
model, where a particle hops from site to site in the pres- 
ence of disorder, and is also allowed to escape the system 
from any site. When the wavelength of the particle is 
comparable with the sample size, an effective long-range 
hopping is created between the sites. This coupling can 
induce the ST, which affects in a non-trivial way the 
transport properties of the system. Similar models of 
quantum transport with coherent dissipation have been 
already considered in the literature [13], but a detailed 
analysis of the interplay of localization and superradi- 
ance has been lacking. This interplay has been recently 
analyzed in Ref. [12], but there the particle was allowed 
to escape only from the end sites, while in the situation 
analyzed in this work, all sites are coupled to the external 
environment. This situation occurs in many important 
physical situations, such as in cold atoms, where a single 
photon is injected in the atomic cloud or in quantum 
dots [15] , 

Intrinsic disorder and opening to the environment have 
opposing effects: while disorder tends to localize the wave 
functions, the opening tends to delocalize them, since it 
induces a long range interaction. The aim of this paper 
is to study the interplay of disorder and opening, and 
the relation to superradiance. We show that while below 
the ST, all states are affected by the disorder and the 
opening in a similar way, above it, the effects are quite 
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different for superradiant and subradiant subspaces, the 
latter being more affected by disorder than the former. 



In Sec.|njwe introduce the model, in Sec. Ill we analyze 
the ST in our system, and in Sec. |IV| we present our main 
numerical results, which we partly justify in Sec. |V| using 
perturbation theory. Finally in Sec. |VI| we present our 
conclusions. 



II. MODEL 

Our starting point is the standard one-dimensional An- 
derson model [HH], for the motion of a particle in a disor- 
dered potential. The Hamiltonian of the Anderson model 
can be written as: 

N N-l 

Ho = '£E j \j)(j\ + nYl (m + il + li + 1)01) . (i) 

3=1 3=1 

where Ej are random variables uniformly distributed in 
[-W/2, + W/2], W is a disorder parameter, and fl is the 
tunneling transition amplitude (in our numerical simu- 
lations we set = 1). For W — the eigenstates are 
extended and we have for the eigenvalues: 
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and the eigenstates: 
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where q = 1, N is a quantum number and j = 1, N 
is a discrete coordinate. In this case, the eigenvalues lie 
in the interval [— 2fl, 20] , so the mean level spacing can 
be estimated as D = 40/iV. For W ^ 0, the eigenstates 
of the one-dimensional Anderson model are exponentially 
localized on the system sites, with exponential tails given 
by 1^0)1 ~ exp(— \j— io|/£); where for weak disorder, the 
localization length £ can be written as: 



£«96 (l-{E/2n) 2 ) 
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For E = 0, Eq. Q has to be modified and we have 
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The phenomenon of Anderson localization was studied 
in a closed disordered chain, while in our case we can vary 
the degree of openness of the system. In particular we 
consider the model in which all the sites are coupled to a 
common channel in the continuum, with equal coupling 
strength 7. This situation can arise when the wavelength 
of the decaying particle is much larger than the size of 
the system. This results in a coherent dissipation, which 
differs from the usual dissipation where every site decays 
independently to a different channel in the continuum. A 
comparison between these two different mechanisms will 
be the subject of a future work. The continuum coupling 



can be taken into account with the aid of an effective 
non-Hermitian Hamiltonian [12 , which in general can be 
written as, 

H eG (E) = H + A(E) - iQ{E) , 

where H is the Hermitian Hamiltonian of the closed 
system decoupled from the environment and A (£7) and 
Q{E) are the induced energy shift and the dissipation, 
respectively. Neglecting the energy dependence and the 
energy shift we have 



{H e g)ij — {Ho)i 
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where A\ are the transition amplitudes from the discrete 
states i to the continuum channels c. 

In the case under study, we have only one decay chan- 
nel, c = 1, and all couplings are equal, so that A\ = ^7. 
Thus the effective Hamiltonian can be written as: 



Ho - i^Q > 



(6) 



where Hq is the Anderson Hamiltonian with diagonal dis- 
order, Eq. |l]), and Qij — 1 Vi,j. 

In order to study the interplay of Anderson localiza- 
tion and superradiance we analyze the participation ratio 
(PR) of the eigenstates of -ff c ff, defined as, 

pr = (e7kW)' (7) 

where the average is over disorder. For completely ex- 
tended states we have PR = N and for completely local- 
ized states we have PR = 1. 



III. SUPERRADIANCE TRANSITION 

ST can be analyzed by studying the complex eigenval- 
ues £ r = E r — iT r /2 of Heff defined in Eq. As the 
coupling between the states and the continuum increases, 
one observes a rearrangement of the widths r r . ST is ex- 
pected to occur for (T)/D ~ 1 [HJ [T^] . The average width 
(r) is 7, so we can define 

k = i/D (8) 

as the parameter controlling the coupling strength to the 
continuum. In the deep localized regime where disorder 
is strong (W 3> 0) we can write D s» W/N, so that the 
effective coupling strength can be written as: 

K = 

W 

In Fig. [T] we show that ST occurs at K ~ 1 for different 
values of W/fl and N. 

For k 3> 1, we can treat the matrix Q as the lead- 
ing term in Eq. (km!, and Ho as a perturbation. The 
superradiant state ^SR) is given to zeroth order by the 
only eigenstate of Q with nonzero eigenvalue: \d) — 
(1, 1) T , and the energy of \SR) is evaluated at first 
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FIG. 1: The average width of the N — 1 subradiant states, 
normalized by the mean level spacing D, versus the effective 
coupling strength k for different values of N and W, and 
= 1. Here we average over 100 disordered configurations. 
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and Ei are the random diagonal elements of Ho. Aver- 
aging over disorder and taking into account that Ei are 
distributed uniformly in [-W/2, W/2] we obtain, 

N -I 

(e)=2ft— — (11) 
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These results agree with our numerical simulations for 
different values of N and allow one to know the position 
in the energy band of the superradiant state in the limit 
k> 1. From Eq. (Ill we deduce that the mean energy 

(e) 



of the superradiant state is independent of W. 



IV. NUMERICAL RESULTS 

In order to study the interplay of superradiance and 
disorder we have analyzed the PR of the eigenstates of 
the non-Hermitian Hamiltonian, Eq. As explained in 
the previous section, as the coupling with the continuum 
is increased we have the formation of one superradiant 
state (the one with the largest width) and N — 1 subra- 
diant ones. In Fig. [2] (upper panel) we analyze the PR 
as a function of k for the states that become subradiant 
for k > 1, and in Fig. [2] (lower panel) we analyze the 
case of the state with the largest width, which becomes 



superradiant for k > 1. As the opening, determined by 
the parameter k, increases, the PR of both superradiant 
and subradiant states increases, showing that the open- 
ing has a delocalizing effect. But the consequences of 
the opening are very different for superradiant and sub- 
radiant states. For the latter, the PR reaches a plateau 
value above the ST (k « 1), which is slightly higher than 
the PR for k <C 1. Moreover on increasing the disorder, 
the PR of the subradiant states decreases, both below 
and above the ST, see Fig. [2] upper panel. The situation 
is different for the superradiant states. Above the ST 
these states become completely delocalized (PR = N) 
and their derealization is not affected by an increase in 
W, see Fig. [2] lower panel. 
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FIG. 2: The participation ratio PR is shown as a function 
of k for different disorder strengths. In the upper panel we 
consider states with —1.5 < E/Q < —0.5, which become sub- 
radiant for large k, while in the lower panel we consider the 
state with the largest width, which corresponds to the super- 
radiant state for large k. Here N = 100, Q — 1, and the PR 
is averaged over 4000 disorder realizations. 

We now look more closely at how the subradiant and 
superradiant states are affected differently by increasing 
the disorder strength W. In Fig. |3j we consider the case 
of N = 100 and 7 = f2 = 1. For small disorder we have 
D w 4n/N, so that 

k = 7/Z) = 1 N/m k 25 > 1 . 

This implies that we are in the superradiant regime. 
Moreover for sufficiently small disorder we have that 
the localization length is larger than the system size, 
£ w 100 £l 2 /W 2 > N, so that both superradiant and 
subradiant states are delocalized. For larger disorder, 
here W > 1, we enter the localized regime, for which 
£ < N. In this regime the PR of the subradiant states 
starts to decrease, while the PR of the superradiant state 
remains unchanged (PR = N), signaling a superradiant 
state that remains completely delocalized. As we increase 
disorder further, k decreases according to Eq. ([9|. The 
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ST occurs at If k 7 AT, here W « 100, and above this 
value the superradiance effect disappears. Summarizing, 
we have a critical value of disorder (W « 100 indicated as 
a full vertical line in Fig. [3]) separating the superradiant 
regime (k > 1), from the non-superradiant one (k < 1). 
Only for W > 100, i.e., below the ST, do the superradi- 
ant states begin to localize, and, for very large disorder, 
corresponding to very small k, they behave the same as 
the subradiant states. 
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FIG. 3: The participation ratio is shown as a function of the 
disorder strength W. Open circles stand for the subradiant 
states, while full circles indicate the superradiant state. Each 
point is obtained by averaging over 100 disorder realizations 
for the superradiant state, while for the subradiant states, 
an additional average over all the subradiant states is per- 
formed. The right and left vertical lines indicate the ST and 
the derealization transition, respectively. Here N = 100 and 

7 = n = 1. 



V. DISCUSSION 

In this Section we will justify (using perturbation the- 
ory) and briefly discuss the interesting results presented 
previously: for small k (below the ST) all the states are 
affected in a similar way by the opening and disorder, 
while for large n (above the ST), the superradiant states 
remain completely delocalized, independent of the degree 
of disorder, while the subradiant states are still sensitive 
to disorder, and their PR decreases with increasing dis- 
order. 



A. Perturbative approach for k< 1 

In the limit k«1, the eigenstates of H c g at first order 
in perturbation theory can be written as: 



1 
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where |n°) are the eigenstates of the closed system, i.e., of 
the Anderson model. Of course, the perturbation expan- 
sion makes sense only when each coefficient in the sum 
in Eq. ( 13 ) is much less than one. This cannot be true 
in general since the eigenvalues E n o are random numbers 
uniformly distributed in the interval [-W/2, W/2}. Thus 
perturbation theory cannot be applied tout court, but 
only for those states whose energies are not too close. 

This simple observation has deep consequences for the 
structure of the eigenstates. Indeed we observe numer- 
ically that on the one hand many single-peaked eigen- 
states become double- or multiple-peaked as 7 increases, 
while on the other hand, they all develop a constant 
plateau proportional to (7/W) 2 , see Fig. [4] 

This last fact can easily be explained using first-order 
perturbation theory as given by Eq. (|13|): in the deep 



localized regime W ^> fl, the matrix elements 
are of order unity and the average distance between two 
random energies is W/3, so the typical coefficients (feo|n) 
in Eq. (13) are ~ 7/W. Furthermore, the mean level 
W/N, and thus the few largest coefficients 
1 N/W ~ k (using Eq. @). 
Thus for weak opening (k <C 1), the typical eigenstate 
consists of a single Anderson model eigenstate with a 
0(k 2 ) admixture of other states, and therefore the typical 
PR for small k differs only by 0(k 2 ) from the PR of the 
Anderson model. 



spacing is D 
in Eq. ( 13 ) are typically 
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FIG. 4: The averaged probability distribution of all eigen- 
states of the non-Hermitian Hamiltonian that are strongly 
peaked in the middle of the chain is shown for different cou- 
pling strength 7 and disorder strength W, as indicated in the 
caption. Specifically, we average over all eigenstates having a 
probability > 0.9 at the site n = N/2 + 1 in order to avoid 
double-peaked states, and also average over disorder. More- 
over, to reduce fluctuations, we average the logarithm of the 
probability distribution. In all cases we fix TV = 100 and 
£1 = 1. Dashed horizontal lines are proportional to (7/W) 2 
in agreement with the perturbative approach. 

As already remarked previously, the perturbative ap- 
proach cannot always work, because for arbitrarily small 
k there is a small but finite probability that two energy 
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states are too close together. This clustering behavior 
has important consequences for the localization proper- 
ties. Specifically, since the nearest-neighbor level spacing 
distribution of uniform random numbers E n o is Poisso- 
nian: P(s) — (1/D) e~ s ' D , where s is the energy differ- 
ence between nearest-neighbor levels and D = W/N is 
the mean level spacing, we can evaluate the probability 



to have two levels closer than 7/2 as 1 — e 



-7/2x1 



K /2for 



small k. This means that there are kN states out of N, 
for which perturbation theory cannot be applied. When 
this happens, the Anderson states mix strongly and the 
PR increases by an O(l) factor. Thus, even though this 
behavior is rare, it makes an O(k) contribution to the 
average PR of the weakly open system, which exceeds 
the 0(k 2 ) contribution from the typical states. Indeed 
the average PR can be evaluated as follow: 



PR 



NkPR 2 + (1 - k)NPRi 
AT 



PR 1 + n(PR 2 -PR 1 ) 



where PR\ and PR 2 refer to the PR of the states for 
which perturbation theory can and cannot be applied. 
Since PR X ~ PR(j = 0) + 0(k 2 ), and PR 2 ~ 0(1), we 
have that PR(-f)—PR(j = 0) ~ k. The numerical results 
in Fig. [5] confirm that the effect of the opening on the PR 
grows as k, instead of the re 2 growth predicted by first- 
order perturbation theory. Here we present the average 
(over disorder) of PR{"f) — PR{l = 0), as a function 
of k = Nj/W for fixed disorder strength and different 
values of the system size. In any case this is quite a 
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FIG. 5: The average increase in the participation ratio, com- 
pared with the closed system, is calculated as a function of k, 
for fixed disorder W = 20 and different system sizes N as indi- 
cated in the legend. In each case the average is performed over 
50000 different eigenstates. The line is PR{^) - PR{0) = 2k. 

delicate point and we postpone its full analysis to a future 
work. 



B. Perturbative approach for k ^> 1 

In the limit k » 1 we consider two cases. First, we 
consider the situation where the nearest neighbor tun- 
neling coupling is 51 = 0, in which case we can follow the 
approach explained in Ref. |17j . This approach will be 



very useful also for the case 7^ 0, which we treat below. 



fi = and k > 1 



If f2 = the Anderson Hamiltonian is diagonal in the 
site basis \j) with eigenvalues Ej distributed uniformly 
in the interval [-W/2, W/2]. The eigenstates of the non- 
Hermitian part —i\Q of the effective Hamiltonian are 
\d) = ^=(1,...,1) T (the superradiant state) with eigen- 
value —i^N, and N — 1 degenerate eigenstates with 
eigenvalue (the subradiant states). We will choose \n) 
in a convenient manner later. Following Ref. [17] we can 
rewrite H e g in the basis of these eigenstates using the 
transformation matrix V, which has as its columns the 
eigenstates of Q: 



H cS = V T H V - i|v T QV : 



-i|N h T 
h H 



(14) 



Here h is a vector of dimension N — 1 with components 
1 N 

-=2 B i0'lM>, (15) 
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while the matrix elements of the (N — 1) x (N — 1) sub- 
matrix H are 

N 



(16) 



Now, we can diagonalize H, 

N 

Hfiv =J2 E i(»\3)(jW) = = e^W). (17) 

i=i 

Following Ref. 17J we obtain 
1 h 
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M \ ^ 



li>, (is) 



where the normalization coefficients are given by 
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In the limit k ^> 1, the eigenstates \fi) of the non- 
Hermitian part of H e g are also eigenstates of H c s . Since 
{d\n) = we have, 



N 

3=1 ^ 



E, 



0. 



(20) 



Therefore each eigenvalue of H lies between two neigh- 
boring levels E n , so that the values are also confined 
in the interval [-W/2, W/2]. 

Let us now estimate the magnitude of the mixing ma- 
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trix elements h^. To do this we compute 



h-h = 



N-l N N 
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and using the completeness relation ^„_j 1 0'|m}(mN) = 
— l/N we have 



h-h = (E 2 ) - (£ 2 ) = AE 2 



This leads to 



AE 



W 



y/W^l ^12(N - 1) 



(22) 



(23) 



Each eigenstate in Eq. (18) is a superposition of all 



the site states \ j) with amplitudes r— M ~ jvfg 
that depend only on the energies Ej and not on the site 
positions j. Nevertheless, each state is quite local- 
ized, since the amplitudes are of order unity for the O(l) 
number of sites whose energy is within a few mean level 
spacings of e (i.e., when |e M — Ej\ ~ D = W/N), and 
small otherwise. This small value of the PR for the sub- 
radiant states should be compared with PR = TV of the 
superradiant states. 

The values obtained above for the subradiant and the 
superradiant states correspond to zeroth-order perturba- 
tion theory. On the other hand first-order perturbation 
theory gives: 
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where are random coefficients with (r 2 ) = 1. We 
see that the exact superradiant state \SR) is a combi- 
nation of the unperturbed superradiant state \d) and a 
small admixture of the unperturbed subradiant states 
|/Lt), and the mixing probability decreases as 1/k 2 for large 
k. Similarly, the admixture of the unperturbed superra- 
diant state \d) in each exact subradiant states \SUB^) 
decreases as 1/(k 2 N). This shows that PR « N for the 
superradiant state and PR ~ 1 for the subradiant states 
when k> 1. 



2. !1 / and k » 1 

As a first step we write the Anderson Hamiltonian Ho 
in terms of its eigenstates \n). Obviously the form of \n) 
will depend on the degree of disorder W. In the following 
we limit our considerations to the large disorder regime, 



so that in the basis of the eigenstates of Hq, the matrix 
elements of Q remain of order one, Q nm ~ 1, and we 
can use the results of Sec. |VB 1[ with the site states and 
energies |j) and Ej replaced by the Anderson eigenstates 
and eigenenergies \n) and E n . 

In Fig. [3] we see that for n > 1 (corresponding to 
W < 100), the superradiant state remains unaffected by 
the increase of disorder, while the subradiant states be- 
come more localized as the disorder strength is increased. 
The results of the previous section can be used to under- 
stand this strongly asymmetric behavior of the PR be- 
tween the subradiant states and the superradiant state. 
Indeed at zeroth order in perturbation theory we can see 
that the superradiant state \SR) \d) is completely delo- 
calized, PR — N, while subradiant states \SUB^) ps \fi) 
become more and more localized as we increase disorder. 



Specifically, the site states \j) in Eq. ( 18 1 are replaced 
with Anderson eigenstates |n), with localization length 
£ oc l/W 2 . This difference persists in first-order pertur- 
bation theory, since the mixing probability between the 
super- and sub-radiant states decreases as 1/k 2 for large 
k, see Eq. pi| . 

Our perturbative approach justifies the results pre- 
sented in Fig. [3j where we can see that the subradiant 
states become increasingly localized as we increase dis- 
order. At the same time Fig. [3] shows that the superra- 
diant state remains completely delocalized even as W is 
increased, until we reach the value W 142.8 (n = 0.7) 
where we find numerically that the ST takes place. The 
perturbative approach shows that superradiant states are 
much less sensitive to disorder because their complex en- 
ergies are at a distance greater than = Wk/2 from 
the subradiant states. 



VI. CONCLUSIONS 

We have studied a 1-D Anderson model with all sites 
coupled to a common decay channel (coherent dissipa- 
tion). Our main motivation was to understand the in- 
terplay of opening and disorder. Increasing the disorder 
tends to localize the states. Increasing the opening, on 
the other hand, reduces the degree of localization, and in 
particular induces a superradiance transition, with the 
formation of a subradiant subspace and a superradiant 
state completely delocalized over all the sites. Our re- 
sults show that, while for small opening all the states 
tend to be similarly affected by the disorder, for large 
opening the superradiant state remains delocalized even 
as the disorder increases, while the subradiant states are 
much more affected by disorder, becoming ever more lo- 
calized as the disorder increases. We have explained these 
effects qualitatively, mainly guided by perturbation the- 
ory. Indeed we have shown that the superradiant state is 
not affected by disorder, since its energy is very distant, 
in the complex plane, from the energies of the subradiant 
states. 
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